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ABSTRACT 

The  MODTRAN5  radiation  transport  (RT)  model  is  a  major  advancement  over  earlier  versions  of  the  MODTRAN 
atmospheric  transmittance  and  radiance  model.  New  model  features  include  (1)  finer  spectral  resolution  via  the  Spectrally 
Enhanced  Resolution  MODTRAN  (SERTRAN)  molecular  band  model,  (2)  a  fully  coupled  treatment  of  auxiliary  molecular 
species,  and  (3)  a  rapid  and  high  fidelity  multiple  scattering  (MS)  option.  The  finer  spectral  resolution  improves  model 
accuracy  especially  in  the  mid-  and  long-wave  infrared;  the  auxiliary  species  will  permit  the  addition  of  any  or  the  entire  suite 
of  HITRAN  species,  along  with  default  and  user-defined  profile  specification;  and  the  MS  option  makes  feasible  the 
calculation  of  Vis-NIR  (Near  InfiraRed)  databases  that  include  high-fidelity  scattered  radiances. 

Keywords:  MODTRAN,  remote  sensing,  hyperspectral,  imaging,  band  model,  infirared,  radiative  transfer,  molecular 
transmittance 

1.  INTRODUCTION 

MODTRAN'’^  has  served  as  the  U.S.  Air  Force  (USAF)  standard  moderate  spectral  resolution  radiative  transport  model 
for  wavelengths  extending  from  the  thermal  InfiraRed  (IR)  through  the  visible  and  into  the  ultraviolet  (0.2  to  10,000.0  pm). 
The  MODTRAN  1  cm''  statistical  band  model  was  developed  collaboratively  by  Spectral  Sciences,  Inc.  and  die  USAF 
Research  Laboratory,  and  it  provides  a  fast  alternative  (100-fold  increase  in  speed)  to  the  USAF  first  principles  and  more 
accurate  line-by-line  (LBL)  radiative  transport  models,  FASCODE^  and  FASCODE  for  the  Environment  FASE^. 
Comparisons  between  MODTRAN  and  FASE  spectral  transmittances  and  radiances  show  agreement  to  within  a  few  percent 
or  better  in  the  thermal  IR.  MODTRAN4  includes  flux  and  atmosphere-scattered  solar  calculations,  essential  conponents  in 
analysis  of  near-IR  and  visible  spectral  region  data  that  are  not  readily  generated  by  LBL  models. 

MODTRAN4  and  its  predecessors  have  been  used  extensively  over  the  last  quarter  century  in  the  design  and  analysis  of 
broadband,  multiband,  and  short-wave  IR/Visible  hyperspectral  imaging  sensors.  However,  conventional  interferometers 
and  many  state-of-the-art  and  next-generation  hyper-  and  ultra-spectral  imaging  sensors  working  in  the  long-  and  mid-wave 
IR  will  operate  at  finer  spectral  resolution  than  MODTRAN4  provides. '  These  sensors  will  also  produce  increased  data 
quality,  quantity  and  utility.  This  has  served  as  a  primary  impetus  for  development  of  MODTRAN5.  The  long-wave 
requirement  for  finer  spectral  resolution  led  to  the  reformulation  of  the  MODTRAN  band  model  equivalent  width  and  line 
tail  absorption  algorithms.  The  result  of  that  work  is  the  0.1  cm''  SERTRAN  (Spectrally  Enhanced  Resolution  MODTRAN) 
molecular  band  model.  The  finer  spectral  resolution  model  provides  high  quality  fine  resolution  predictions.  Development 
of  SERTRAN  has  also  led  to  improved  coarser  resolution  statistical  band  models  by  providing  a  source  for  tuning  the  coarse 
spectral  resolution  databases  and  by  refining  the  transmittance  algorithms  within  MODTRAN. 

The  introduction  of  the  auxiliary  species  option  enables  MODTRAN  users  to  simulate  the  impact  of  trace  molecular 
gases.  User  input  of  both  profile  and  optical  data  is  required.  Auxiliary  database  files  for  all  the  HITRAN  molecular  species 
not  already  included  in  MODTRAN  are  being  generated.  In  addition,  the  default  vertical  profiles  for  these  HITRAN-specific 
gases  will  be  made  available.  It  is  anticipated  that  a  user-defined  option  for  other  gases  and  their  profiles  will  eventually 
come  on  line. 
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A  new  multiple  scattering  option  automatically  performs  DISORT^  multiple  scattering  calculations  at  a  few  fixed 
spectral  points  and  then  spectrally  scales  fast  ISAACS®  2-stream  results.  In  the  visible  and  NIR  transparent  window  regions, 
where  multiple  scattering  is  the  dominant  radiation  soiuce,  the  scaling  dramatically  improves  the  accuracy  of  predicted 
radianceSi  TTiese  MODTRAN5  advances  are  described  in  the  following  sections. 

2.  THE  SERTRAN  BAND  MODEL 


Narrowing  the  band  model  spectral  resolution  changes  the  fundamental  character  of  the  band  model.  The  average  half¬ 
width  of  molecular  transitions  near  sea-level  is  of  order  0.07  cm'',  which  is  comparable  to  the  SERTRAN  spectral  bin  width. 

As  illustrated  in  Figure  1,  the  1.0  cm''  statistical  band  model  calculates  the  absorption  of  molecular  lines  whose  line  center 
regions  lie  almost  entirely  within  the  spectral  bin.  A  much  larger  fraction  of  each  molecular  line  falls  outside  of  the  spectral 
bin  at  the  finer  spectral  resolution  of  SERTRAN.  Reformulation  of  the  molecular  band  model  has  therefore  required 
improved  treatment  of  both  line  center  and  line  tail  absorption.  Line  tail  absorption  must  be  modeled  near  line  centers,  and 
the  finite-bin  single-line  equivalent  width  used  to  calculate  line  center  absoiption  is  not  simply  a  small  perturbation  of  the 
total  single  line  equivalent  width. 


Figure  1.  Molecular  transition  spectral  absorption  near  the  ground:  (a)  Spectral  features  within  a  MODTRAN  1.0  cm''  bin;  (b)  Spectral 
features  within  a  SERTRAN  0. 1  cm''  bin. 


The  band  model  line  center  absorption  within  a  spectral  bin  is 
defined  as  the  in-band  absorption  from  all  molecular  transitions 
centered  in  that  bin.  Figure  2.  LBL  models  calculate  this  in-band 
absorption  by  oqdicitly  determining  the  spectral  absorption  of  each 
molecular  line  on  a  fine  spectral  grid  and  then  integrating  the  resulting 
spectra.  In  a  band  model  approach,  the  in-band  absorption  is 
approximated  by  a  statistical  representation  of  the  number  of  lines  and 
their  strength,  location  and  overlap.  Temperature-dependent  band 
model  data  are  computed  from  a  molecular  transition  line  atlas  such  as 
HITRAN’.  These  parameters  define  an  effective  single-line 
absorption  strength  for  the  interval  and  an  effective  number  of  lines. 

The  finite  spectral  bin  single-line  Voigt  equivalent  width  is 
computed  to  determine  the  absorption  of  the  effective  average  line.  It 
is  calculated  as  the  difference  between  the  total  equivalent  width  and 
the  absorption  due  to  the  line  tails  falling  outside  of  the  spectral  band. 
The  line  tail  subtractions  are  performed  for  an  off-centered  line  with 
the  distance  from  spectral  bin  edge  determined  from  the  weak-line 
Lorentz  limit.  Offsetting  the  location  of  the  effective  line  from  the  bin 
center  lowers  its  in-band  absorption,  thereby  producing  a  more 
representative  result  for  a  line  randomly  located  in  the  bin. 
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Figure  2.  Comparison  of  Band  Model  and  Line-By- 
Line  Approaches.  LBL  models  calculate  high  spectral 
resolution  absorption  (solid  curve).  Band  models 
statistically  determine  the  average  integrated  strength 
(dash  line)  based  on  line  strength,  line  density,  and  line 
width  parameters. 


Line  tails  are  not  simply  inversely  proportional  to  square  of  the  line  center  displacement  (the  MODTRAN4  assumption) 
when  the  spectral  bin  width  (Av)  is  0.1  cm"'  and  the  effective  average  line  is  positioned  Av/4  from  the  bin  edge.  Neither  are 
the  line  tails  simply  Lorentzian.  In  SERTRAN,  the  integrated  Voigt  line  tails  are  computed  from  a  newly  derived  exact 
expansion: 


W,  =  Su 


Cl) 


Here,  A  is  the  distance  from  line  center,)^  is  the  normalized  Voigt  line  shape  function,  Su  is  the  line  strength  and  molecular 
column  product,  /„  are  modified  Bessel  fimctions  of  integer  order,  and  flie  V„  are  Fourier  expansion  coefficients: 


K 


Details  of  the  derivation  of  these  expressions  have  been  provided  elsewhere* 


MODTRAN  band  model  line  tail  absorption  for  a  spectral  bin  is  defined  as  the  absorption  from  molecular  transitions 
centered  outside  of  the  bin  but  no  more  than  25  cm  '  from  band  center.  Contributions  from  beyond  25  cm'  are  modeled  as 
continua.  In  MODTRAN5,  die  line  tail  spectral  absorption  is  fit  to  tenperature  and  pressure  dependent  Pad6  approximants. 
Five  parameters,  ko,  *i,  ki,  Xj  and  Xz,  are  used  in  the  Fade  fits: 

_  kfi  +  k^  Sy  +  ^2  .  g  _  ^  ^cat  ^  (3) 

P//»  ~  l  +  x^S^  +  x^S;  V  "  AW2 

In  this  parameterization,  K:*,  is the  central  frequency  of  the  spectral  bin  and  the  range  of  the  normalized  spectral  variable  Sy  is 
from  -1  to  +1.  Since  die  Lorentz  line  shape  to  leading  order  is  proportional  to  pressure,  die  absorption  coefficient  over 
pressure  P  normalized  to  1  atm  (Fo)  is  fit.  This  form  for  the  tail  contributions  is  exact  in  the  limit  of  a  single  Lorentz  tail. 
The  five  Fade  parameters  are  determined  from  the  summed  line  tail  spectral  absorption  coefficients  at  vten  and  edges  of  the 
spectral  band,  the  spectral  frequency  derivative  of  the  absorption  coefficient  at  and  die  integral  of  the  spectral  absorption 
coefficient  over  the  spectral  band. 

Figure  3  illustrates  an  exair^le  of  a  Fade  spectral  fit  to  line  tail  cpiitnbutions.  In  this  case,  hundreds  of  molecular 
transitions  contribute  to  the  305K  and  1  atm  pressure  HjO  line  tail  absorption  coefficient  spectrum  between  200.0  and  200.1 
cm'',  the  thick  dashed  curve.  As  is  generally  tme,  the  spectral  fit  lies  directly  on  top  of  the  high-resolution  spectrum 

The  Fade  p^ameter  database  is  generated  for  6  evenly  spaced  temperatures  from  180  to  305K  and  for  pressures  of  1.0 
and  0.1  atm.  In  MODTRAN,  a  parameter  defines  the  number  of  monochromatic  line  tail  absorption  coefficient  calculations 
per  spectral  bin;  it  is  nominally  set  to  5.  The  line  tail  optical  depth  is  modeled  as  varying  linearly  between  spectral  points  for 
computing  transmittances.  The  k/P  interpolations  are  linear  in  ten:5)erature  and  in  pressure  squared. 

MODTRAN5  spectral  transmittances  from  H2O  absorption,  nadir-viewing  through  the  relatively  wet  Mid-Latitude 
Summer  (MLS)  model  atmosphere,  were  conqjared  to  calculations  from  the  Line-By-Line  Radiation  Transport  Model 
(LBLRTM).  A  common  database  of  molecular  line  data,  the  HITRAN2k’* with  2001  updates,  was  used  for  these 
calculations.  The  spectral  range  extended  from  0.88  to  1.22  pm,  spanning  both  the  0.94  and  1.38  pm  H2O  bands.  The 
results,  illustrated  for  a  small  spectral  subset  in  Figure  4,  are  typical.  The  overlay  of  MODTRAN5  md  LBLRTM  0.2  cm' 
resolution  spectra  show  excellent  agreement.  The  more  detailed  residual  curves  (offset  by  1.10  in  Figure  4)  show  that 
deviations  as  large  as  0.07  can  arise.  These  substantial  deviations  are  a  result  of  the  band  model  formulation,  which  requires 
that  line  centers  be  shifted  by  as  much  as  0.025  cm  '  from  the  edge  of  the  0.1  cm''  spectral  bands.  The  line  center  shifts  in 
turn  produce  a  shift  of  absorption  between  neighboring  bands.  As  the  results  are  further  degraded,  residual  decrease  and 
there  is  generally  little  or  no  biasing  to  too  much  or  too  little  integrated  absorption. 


Frequency  (OTi^) 

Figure  3.  High  Spectral  Resolution  Line  Tail  Absorption  Curve 
(thick  dash)  and  Pad6  Approximant  Fit  (thin  line)  for  Lorentzian 
HjO  Lines  at  305K  and  latm  Pressure.  The  spectral  curves  from 
largest  contributors  are  also  illustrated  (lines  with  diamond  for 
lower  frequency  transitions  and  lines  with  squares  for  higher 
frequency  transitions). 
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Figure  4.  Intercomparison  of  0.2  cm"'  spectral  transmittances 
for  H2O  between  MODTRAN5  and  LBLRTM.  An  offset  by  1 . 1 
has  been  added  to  the  MODTRAN5  -  LBLRTM  residuals. 


3.  AUXILIARY  SPECIES 


In  MODTRAN5,  auxiliary  species  molecular  densities  can  be  defined  via  user-specified  atmospheric  vertical  profiles  or 
with  the  ‘default’  profiles  currently  provided  with  FASCODE/FASE.  For  each  auxiliary  species,  MODTRAN5  accepts  input 
of  either  MODTRAN-like  molecular  band  model  pmmeters  or  absorption  cross-sections,  similar  to  those  aheady  en:q)loyed 
for  the  O3  Chappuis  bands  (non-temperature  dependent)  or  CFCs.  The  band  model  provides  greater  accuracy  if  optical 
depths  are  substantial,  but  absorption  data  is  often  more  readily  available  for  non-HITRAN  species.  Any  of  the  NIST/EPA 
Gas  Phase  Infrared  Database  JCAMP-DX"  (Joint  Comriiittee  for  Atomic  and  Molecular  Physical  Data  “Data  Exchange”) 
measured  absorption  spectra  as  well  as  other  gaseous  database  spectra  can  be  introduced  However,  because  the  actual 
database  is  the  property  of  the  measurement  laboratories,  MODTRAN5  will  not  be  distributed  with  species  spectral  libraries, 
but  only  the  conversion  utility^. 

As  the  column  amoimt  of  an  auxiliary  species  increases,  the  growth  of  in-band  optical  depth  slows  and  a  band  model  is 
required  to  accurately  compute  absorption.  As  noted,  band  model  data  files  are  being  generated  for  HITRAN  species  not 
already  included  in  MODTRAN.  Alternatively,  band  model  data  can  be  computed  from  either  extensive  compilations  of 
transition  frequencies,  line  strengths,  lower  state  energies  and  collision  broadened  half- widths,  or  from  absorption  spectra  as 
a  function  of  molecular  column  density  for  a  series  of  temperatures. 

The  mechanism  for  introducing  auxiliary  species  is  similar  to  that  for  MODTRAN  standard  atmospheric  molecular 
species  when  user-defined  profiles  or  radiosonde  data  are  being  entered.  For  each  atmospheric  level,  the  density  of  each 
species  is  input  directly  or  is  linked  to  a  default  profile.  Unlike  the  standard  molecules,  the  munber  of  auxiliary  species  is 
read  from  the  profile  table  header.  This  number  is  the  total  number  of  auxiliary  species  in  the  user’s  master  list.  The  input  to 
MODTRAN5  accepts  a  negative  sign  prefix  for  any  species  name  to  indicate  the  species  is  to  be  excluded.  This  input 
format  was  instituted  so  that  flie  input  stream  would  not  have  to  be  completely  reconfigured  each  time  a  species  is  to  be 
added  or  removed.  A  star  (*)  at  the  end  of  species  name  indicates  band  model  rather  than  absorption  coefficient  data  is  to  be 
read.  The  optical  data  for  each  species  is  stored  in  species-dependent  auxiliary  input  files  (SF6.BM  for  a  SFg  absorption 
coefficient  file,  SF6.CBM  for  a  SFs  line  center  band  model  data  file,  and  SF6.TBM  for  a  SFe  line  tail  band  model  data  file). 


To  illustrate  the  auxiliary  species  option,  MODTRAN  band  model  data  files  have  been  generated  for  HCl  diatonnic 
molepule^  specifically,  files  HCICBM  and  HCITBM.  A  simple  profile  containing  a  surface  layer  of  HCl  was  folded  in  the 
MODTRAN  MLS  model  atmosphere  by  setting  the  HCl  density  to  a  constant  value  between  0.00  and  0.99km  and  having  it 
ramp  down  to  zero  over  the  next  20  meters.  MODTRAN  nadir-viewing  spectral  transmittances  and  radiances  were  generated 
for  HCl  concentrations  of  0.1,  1.0  and  10.0  ppm.  For  a  Top-Of-Atmosphere  sensor,  HCl  transmittance  spectra  are  illustrated 
in  Figure  5.  Oh  the  left,  the  envelope  of  R  and  P  branch  transitions  from  the  (0-1)  vibrational  transition  are  displayed.  At  0.2 
cm'*  spectral  resolution,  the  splitting  of  the  HCl^^  and  HCl^^  isotopes  is  resolved.  This  is  illustrated  for  R2  and  P3  transitions 
on  the  right  in  Figure  5. 
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Figure  5.  HCl  Spectral  Transmittances  as  Predicted  by  the  MODTRAN5  using  the  Auxiliary  Species  Option 

The  transmittance  spectrum  of  the  R2  and  P3  vibrational/rotational  lines  for  the  HCl  fundamental  arising  from  the  HCl 
surface  layer,  when  embedded  in  the  MODTRAN  MLS  atmosphere  and  viewed  from  space.  Figure  6,  is  considerably  more 
con^lex  than  the  spectmm  of  the  isolated  species  because  of  spectral  overlap  with  other  molecular  species.  Both  CH4  and 
H2O  have  strong  absorption  lines  in  this  spectral  region.  The  profiles  of  these  con^eting  absorbers  must  be  well 
characterized  if  HCl  is  to  be  correctly  included  in  the  simulation.  There  is,  of  course,  a  feedback  between  the  temperature 
and  the  distribution  of  all  gases  in  a  radiative  transfer  calculation,  whether  modeled  by  a  band  model  like  MODTRAN  or  an 
LBL  code.  In  Figure  6  the  concentration  of  HCl  increases,  while  the  observed  thermal  emission  increases.  If,  however,  the 
surface  temperature  was  not  lowered,  the  effect  of  the  added  gas  would  have  been  to  lessen  the  observed  radiance.  In  effect, 
the  HCl  emission  would  be  occurring  at  a  lower  atmospheric  tenqjerature  than  the  background  surface.  Thus,  the  result  from 
introducing  auxiliary  species  (or  altering  the  concentration  and/or  temperature  of  any  of  the  extant  species)  is  a  complex 
function  dependent  on  contrast  between  surface  and  gas  temperatures. 


Figure  6.  Spectral  Transmittances  and  Thermal  Emission  as  Predicted  by  the  MODTRAN5  for  Nadir-Viewing  through  a  MLS  atmosphere 
with  a  1  km  High  HCl  layer  at  the  Surface.  The  blackbody  ground  temperature  is  lOK  below  the  surface  air  temperature. 


4.  RAPID  ACCURATE  MULTIPLE  SCATTERING 

Previous  versions  of  MODTRAN  provide  twb  algorithms  for  confuting  the  multiple  scattering  conponents  of  the  total 
path  radiance:  an  approximate  ISAACS  2-streaiti  method  and  the  accurate  but  time-consunaing  multi-stream  DISOR.T 
(Discrete  Ordinate  Radiance  Transfer)  method.  DISORT  runs  100  to  1000  times  slower  than  the  2-stream  method.  For 
MODTRAN5,  an  option  has  been  developed  to  spectrally  scale  the  2-stream  calculations  using  DISORT  results  at  a  small 
number  of  wavelengths.  This  scaling  recovers  almost  full  DISORT  accuracy  with  less  than  a  factor  of  two  time  penalty  as 
compared  to  ISAACS. 

The  DISORT  scaling  works  as  follows.  Twelve  spectral  wavelengths  between  2.14  and  0.20  pm  have  been  selected 
lying  outside  of  the  stronger  MODTRAN  atmospheric  molecular  bands,  Figure  7.  These  values  are  hardwired  into 
MODTRAN5.  Molecular  bands  are  avoided  because  multiple  scatter  radiance  contributions  are  less  dominant  in  the 
presence  of  molecular  absorption.  As  an  automated  preprocess  within  MODTRAN,  DISORT  and  ISAACS  inultiple 
scattering  path  radiance  calculations  with  the  user-defined  solar  and  path  geometry  are  performed  at  the  DISpRT  scaling 
wavelengths  encon^assing  die  input  spectral  range.  Two  spectral  correction  curves  are  generated  from  these  calculations: 
one  for  the  path  scattered  radiance  (excluding  single  scatter  solar)  and  the  second  for  flie  downward  diffuse  flux  at  the 
ground.  The  corrections  are  defined  as  the  ratio  of  flie  DISORT  to  ISAACS  values  and  linearly  interpolation  between 
wavelengths.  These  spectral  corrections  are  then  applied  to  the  ISAACS  results  computed  at  the  band  model  spectral 
resolution. 

Figure  8  illustrates  the  benefit  of  scaling.  In  this  figure,  spectral  radiances  are  simulated  for  a  nadir  viewing  fi-om  space 
and  a  40“  sun.  The  atmosphere  was  set  to  the  MODTRAN  mid-latitude  summer  with  its  relatively  high  water  column. 
Aerosol  loading  was  also  relatively  high  with  the  surface  meteorological  range  (visibility)  input  of  23-km  and  Navy  Mafrtime 
aerosol  type.  The  surface  albedo  is  zero,  so  only  the  path  niultiple  scattered  radiance  is  scaled.  The  spectral  radiance 
predictions  in  the  visible  spectral  region  (400-700hm)  with  ISAACS  multiple  scattering  are  —10%  below  die  resiilts  with  8- 
stream  DISORT.  inh-ndiiiring  the  scaling  algorithm  essentially  eliminates  these  residuals.  In  the  Near-IR,  the  overall 
performance  of  scaling  is  mixed;  Outside  of  the  Strong  molecular  bands  the  sealirig  significantly  hmiroves  the  radiaiiCe 
results,  but  within  molecular  bands  die  scaling  produces  too  much  radiance.  In  these  rdgioris,  best  results  are  Still  obtained  by 
performing  the  rigorous  DISORT  calculations. 


Wavelength  (pm)  Wavelength  (nm) 

Figure  7.  Spectral  Grid  for  DISORT  Scaling  of  ISAACS  2-  Figure  8.  Results  for  DISORT  Spectral  Scaling  of  ISAACS 

Stream  Multiple  Scattering.  Wavelengths  are  selected  outside  of  Radiances  with  a  Black  Surface  (Albedo  =  0). 

Strong  atmospheric  molecular  bands. 


In  a  second  example,  the  results  of  scaling  are  illustrated  with  a  fully  reflective  groimd  sxuface.  Figure  9.  As  expected, 
the  spectral  radiances  are  almost  an  order  of  magnitude  higher  than  in  Figure  8.  The  results  of  scaling  remain  good  outside  of 
molecular  bands.  Within  these  bands  the  residuals  are  substantial;  for  best  results,  DISORT  should  be  used  within  these 
bands.  The  main  benefit  of  scaling  is  in  analysis  and  retrieval  of  aerosols.  Traditionally,  the  spectral  window  regions  hawe 
been  exploited  to  extract  aerosol  loading  and  type.  This  analysis  is  greatly  improved  with  accinate  multiple  scattering,  which 
the  DISORT  scaling  option  provides. 


Nadir  from  Space  - 1.0  Surface  Albedo 


Figure  9.  Results  for  DISORT  Spectral  Scaling  of  ISAACS  R^iances  with  a  White  Surface  (Albedo  =  1). 

5.  CONCLUDING  REMARKS  AND  ACKNOWLEDGEMENTS 

This  paper  has  demonstrated  that  MODTRAN5  provides  capabilities  not  previously  available  to  the  research  community, 
which  are  expected  to  add  to  the  flexibility  of  atmospheric  radiative  transfer  modeling  for  the  community  at  large.  The 
authors  wish  to  acknowledge  the  generous  support  of  the  Air  Force  Research  Laboratory  (AFRL),  the  National  Imagery  and 
Mapping  Agency  (NIMA),  the  Spectral  Information  Technology  Applications  Center  (SITAC),  the  Department  of  Energy 
(DOE),  National  Aeronautics  &  Space  Administration  (NASA).  Also  thanks  to  Dr.  H.E.  Snell  (Atmospheric  and 
Environmental  Research,  Inc.)  who  provided  line-by-line  calculations. 

REFERENCES 

1.  A.  Berk,  G.P.  Anderson,  P.K.  Acharya,  L.S.  Bernstein,  J.H.  Chetwynd,  M.W.  Matthew,  E.P.  Shettle  and  S.M.  Adler- 
Golden,  “MODTRAN4  User's  Manual,”  Air  Force  Research  Laboratory  Report,  June  1999. 

2.  A.  Berk,  L.S.  Bernstein,  G.P.  Anderson,  P.K.  Acharya,  D.C.  Robertson,  J.H.  Chetwynd  and  S.M.  Adler-Golden, 
“MODTRAN  Cloud  and  Multiple  Scattering  Upgrades  with  Application  to  AVIRIS,”  Remote  Sens.  Environ.  65,  pp. 
367-375, 1998. 

3.  S.A.  Clough,  F.X.  Kneizys,  G.P.  Anderson,  E.P.  Shettle,  J.H.  Chetwynd,  L.W.  Abreu,  and  L.A.  Hall,  “FASCODE3 
Spectral  Simulation,”  Proceedings  of  the  International  Radiation  Symposium,  Lenoble  and  Geleyn,  Deepak  Publishing, 
1988. 

4.  H.E.  Snell,  J.-L.  Moncet,  G.P.  Anderson,  J.H.  Chetwynd,  S.  Miller,  and  J.  Wang,  “FASCODE  for  the  Environment 
(EASE)”,  Proceedings  of  Atmospheric  Propagation  and  Remote  Sensing  IV,  SPIE,  2471,  pp.  88-95, 1995. 

5.  K.  Stamnes,  S.C.  Tsay,  W.J.  Wiscombe  and  K.  Jayaweera,  “Numerically  Stable  Algorithm  for  Discrete-Ordinate-Method 
Radiative  Transfer  in  Multiple  Scattering  and  Emitting  Layered  Media”,  Applied  Optics,  27,  pp.  2502-2509,  1988. 

6.  R.G.  Isaacs,  W.C.  Wang,  R.D.  Worsham  and  S.  Goldenberg,  “Multiple  Scattering  LOWTRAN  and  FASCODE  Models,” 
Applied  Optics  26,'pp.\212-\2i\,\9Zl. 

7.  L.S.  Rothman,  C.P.  Rinsland,  A.  Goldman,  S.T.  Massie,  D.P.  Edwards,  J.-M.  Flaud,  A.  Perrin,  V.  Dana,  J.-.Y.  Mandin, 
J.  Schroeder,  A.  McCann,  R.R.  Gamache,  R.B.  Wattson,  K.  YoShino,  K.  Chance,  K.W.  Jucks,  L.R.  Brown,  V. 
Nemtchinov,  and  P.  Varanasi,  The  HITRAN  Molecular  Spectroscopic  Database  and  HAWKS  (HITRAN  Atmospheric 
Workstation):  1996  Edition,  J.  Quant.  Spectrosc.  Radiat.  Transfer,  60,  pp.  665-710  (1998). 


8.  A.  Berk,  P.K.  Acharya,  L.S.  Bernstein,  G.P.  Anderson,  P.E.  Lewis,  J.H.  Chetwynd  and  M.L.  Hoke,  “Band  Model 
Method  for  Modeling  Atmospheric  Propagation  at  Arbitrarily  Fine  Spectral  Resolution”,  Application  for  United  States 
Letters  Patent,  submitted  April  18,  2001. 

9.  http://atmos.umd.edu/~bobe/LBLRTM/. 

10.  http://www.hitran.com/. 

11.  http://www.icamp.org/. 


